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ABSTRACT 

We present the analysis of simultaneous multi-frequency Very Large Array (VLA) 
observations of 57 out of 61 sources from the "faint" high frequency peaker (HFP) 
sample carried out in various epochs. Sloan Digital Sky Survey SDSS data have been 
used to identify the optical counterpart of each radio source. From the analysis of 
the multi-epoch spectra we find that 24 sources do not show evidence of spectral 
variability, while 12 objects do not possess a peaked spectrum anymore at least in 
one of the observing epochs. Among the remaining 21 sources showing some degree of 
variability, we find that in 8 objects the spectral properties change consistently with 
the expectation for a radio source undergoing adiabatic expansion. The comparison 
between the variability and the optical identification suggests that the majority of 
radio sources hosted in galaxies likely represent the young radio source population, 
whereas the majority of those associated with quasars are part of a different population 
similar to fiat-spectrum objects, which possess peaked spectra during short intervals 
of their life, as found in other samples of high-frequency peaking objects. The analysis 
of the optical images from the SDSS points out the presence of companions around 6 
HFP hosted in galaxies, suggesting that young radio sources resides in groups. 

Key words: galaxies: active - galaxies: evolution - radio continuum: general - quasars: 
general 



> 1 INTRODUCTION 
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Our knowledge of the first stages of the evolution of pow- 
erful radio sources is based on the study of the population 
of high frequency peaking radio sources. In the framework 
of models explaining the evolution of individual radio 
sources, the spectral peak of young radio sources occurs at 
high frequencies. Given their small size, in these sources 
the synchrotron self-absorption (SSA) is a very efi'ective 
mechanism. As the source grows, the peak frequency is 
expected to shift towards lower frequencies as a conse- 
quence of adiabatic expansion. An alternative explanation 
suggests that the spectral peak is due to free-free absorption 
from a ionized med ium enshrouding the radio emission 
l|Bicknell et al.lll997l ). Both scenari os are supported by the 
empirical anti-correlation found bv lO'Dea &: Baum I (|l997l ) 
from the study of samples of compact steep spectrum (CSS) 
and gigahertz-peaked spectrum (GPS) radio sources. The 
former have peak frequencies around a few hundred MHz, 
typical sizes of a few kpc and ages of 10^ - 10^ years, whereas 
the latter have spectra peaking around 1 GHz, typical sizes 
of about 1 kpc or less and ages of 10^ - 10'* years. However, 
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it is worth noting that the consistency between the source 
size and the spectral peak often found in the most compact 
sources s trongly support the synchro t ron self-absorption 
scenario lOrienti fc Dallacasa I l2008al : iTingav fc de Kool I 

120031) . 

High frequency peakers (HFP: [Pallacasa et al.|[2000l ). with 
a spectral peak occurring at frequencies above a few GHz, 
are thus the best candidates to be newly born radio sources, 
with ages between lO'^ - 10"^ years. 

The radio properties of HFPs have been derived by the 
analy si s of the " b right" HFP sample (jPallacasa et al.| 
2000l: iTinti et al.l [2OO5I: lOrienti et all 2006a, 20o3; 
Orienti fc Dallacasa I r2008bl ). In particular, from the 



multi-epoch analysis of their radio spectra it has been found 
that the sample is composed of two different populations. 
One population consists of radio sources that maintain the 
convex spectrum without showing variability, whereas the 
other comprises radio sources that change their spectral 
shape, becoming also flat-spectrum objects, and possessing 
substantial flux density variability. The different spectral 
properties shown by the two populations suggest that the 
former represent young radio sources still in an early stage 
of their evolution, while the latter are beamed objects. 
The analysis of the multi-epoch spectral behaviour has 
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Table 1. VLA observations and configurations 



Date 


Conf 


Proj.ID 


Code 


Sep. 2003 


AnB 


AD488 


a 


Jan.2004 


BnC 


AD494 


b 


Nov.2006 


C 


AO210 


c 


Apr. 2007 


D 


AO210 


d 



proved to be a p owerful tool to discrim i nate between the 
two populations (iTorniainen et al] l2005l : iTinti et al.l l2005l: 
lOrienti et al]|200'rf r In fact, beamed radio sources, although 
usually characterized by a flat and variable spectrum, 
may be selected in samples of high-frequency peaking 
objects when their emission is dominated by a flaring 
knot in a jet. On the other hand, young radio sources are 
known t o be the least variable class of extragalactic radio 
sources (jO'Deal [ 19981 ). However, it must be mentioned 
that in the youngest objects, substantial variability in 
the optically-thick part of the spectrum is expected as a 
consequence of either the source growt h / evolution (e.g. 
J1459-I-3337, lOrienti fc Daflacasa I l2008d ). or changes in 
the possible absorber screen, or a combination of both 
jXingav fc de Kool 2003). 

In this paper we present a multi-epoch analysis based 
on simultaneous multi-frequency VLA data of the radio 
spectra of 57 high frequency peaker s from the "faint" 
HFP sample (|Stanghellini et alj |2009|). This sample was 
selected as the "bright" HFP sample ( Dallacasa et aklliood ) 
by cross-correlating the 87GB survey at 4.9 GHz with 
the NVSS at 1.4 GHz, and including only sources fainter 
than 300 mjy at 4.9 GHz within a restricted area around 
the northern galactic cap (for details on the sample see 
Stanghellini et al. 2009) . The study of the radio properties 
of a sample of faint HFFs is the first step in understanding 
the first stages of the radio source evolution. So far, spectral 
studies have been carried out for the bright HFP sources 
only, and an extension to fainter objects is necessary. The 
evolution models developed so far (e.g . Fanti et all 1 19951 : 



ISnellen et al]l2000l : iKaiser fc Alexandejll997l 'l predict that 
in the earliest stage the radio luminosity progressively 
increases, implying that the youngest objects are likely 
to be found among faint sources. Furthermore, in faint 
HFPs, boosting effects should be less relevant, making the 
contamination from blazars less severe than what found in 
samples of brighter objects. 

Throughout this paper we assume the following cosmol- 
ogy: Ho = 71km s"^ Mpc"\ = 0.27, and Q,k = 0.73 in 
a flat Universe. The spectral index is defined as S{v) oc p^" . 



2 VLA OBSERVATIONS AND DATA 
REDUCTION 

Simultaneous mult i- frequency VLA observations of 57 
out of the 61 sources from the "faint" HFP sample 



IStanghellini et al.|[2009l) were carried out during different 
runs between September 2003 and April 2007 (Table [l|. 
Observations were performed in L band (with the two IFs 
centered at 1.415 and 1.665 GHz), C band (with the two 
IFs centered at 4.565 and 4.935 GHz), X band (with the 
two IFs centered at 8.085 and 8.465 GHz), U band (14.940 
GHz), K band (22.460 GHz), and in Q band (43.340 GHz). 
At each frequency, the target sources were observed for 
about 1 minute, cycling through frequencies. During each 
run, the primary fiux density calibrator either 3C 286 or 
3C 48 was observed for about 3 minutes at each frequency. 
Secondary calibrators were chosen on the basis of their 
distance from the targets in order to minimize the telescope 
slewing time, and they were observed for 1.5 min at each 
frequency, every 20 min. 

The data reduction was carried out following the standard 
procedures for the VLA implemented in the NRAO AIPS 
package. The flux density at each frequency was mea- 
sured on the flnal image produced after a few phase-only 
self-calibration iterations. In the L band it was generally 
necessary to image a few confusing sources falling within the 
primary beam. All the target sources appeared unresolved 
at any frequency. During the observations of a few sources, 
strong RFI at 1.420 and 1.665 GHz was present, and in 
those cases the measurements of the flux density were not 
possible. 

Uncertainties on the determination of the absolute fiux 
density scale are dominated by amplitude errors. Based on 
the variations of antenna complex gains during the various 
observations, we can conservatively estimate an uncertainty 
of ~3% in L, C, and X bands, ^^5% in U band, and ~10% 
in K and Q bands. The rms noise level on the image plane 
is about 0.1 mjy/beam in C and X bands, and about 
0.2-0.4 mJy/beam in L, U, and K bands. In the Q band 
it accounts for 0.4-0.5 mJy/beam, becoming comparable 
to the amplitude calibration errors for sources fainter than 
~20 mJy. Results are presented in Section 4. 



3 OPTICAL DATA 

To determine the optical properties of the sources in 
the faint HFP sample, we complemented the information 
availa ble in the literature w ith that provided by the SDSS 
DR7 (|Abazaiian et al]|2009l '). The optical properties of each 
object (like source extension and magnitude) have been 
carefully inspected beyond the automated procedures in 
the SDSS, in order to unambiguously identify the host (i.e. 
quasar, galaxy, or empty field) of each radio source. 
Of the 57 sources considered in this paper, 12 are identified 
with galaxies with redshift between 0.03 and 0.6; 33 are 
quasars with a higher redshift, typically in the range from 
0.6 to 3.0, while 12 sources still lack an optical counterpart 
(labelled as empty field (EF) in Table [3|. 
Images and optical information have been retrieved by 
means of the SDSS DR7 Finding Chart Tool. In Table E] we 
report the R magnitude, converted in the Johnson-Kron- 
Cousins BVRI system, and the spectroscopic redshift when 
available. An "f" indicates a photometric redshift. 
In Fig. [T]we present th e optical images from the SDSS DR7 
(|Abazaiian et al]|2009l ) of the sources hosted in galaxies. 
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Figure 1. Optical images from the SDSS DR7 of the 12 HFP radio sources identified with a galaxy. On the images we report the source 
name, the redshift (an "f" indicates a photometric redshift). If the source forms a group, we report the redshift of the companion galaxies, 
when available. The field width has been chosen to show a region of about 250 kpc around the galaxy hosting the HFP radio source. 



For the 3 galaxies with available spectroscopic data we show 
the spectrum in Fig. [2] and we summarize the information 
on the main lines in Table (2] 



3.1 Optical properties 

A characteristic arising from Fig. [T] is the presence of 
companions within a projected distance of about 150 - 
200 kpc around 6 HFP galaxies. Although in J0804+5431, 
J1058+3353, J1109+3831, and J1218+2828 this evidence 
comes from photometric information only, in the case of 
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Figure 1. Continued. 



J1530+2705 and J1602+2646 the association is confirmed 
by spectroscopic redshifts, supporting the idea that young 
radio sources reside in groups, as found in other works 
llOrienti et al. I l2006bl : ISnelien et al. l l2002l : IStanghellini et al.l 
liggd ). Therelativelv small redshift of J1530+2705 aUows 
us to identify the spiral morphology of the brightest 
companions where also a bar is clearly visible. 
The galaxies hosting the HFPs are usually the brightest 
elliptical at the group centre. In the case of J0804+5431 
the radio source is hosted by an elliptical galaxy that is 
at a projected distance of about 160 kpc to the north-east 
from the brightest galaxy at the centre of the group and it 
is at a projected distance of about 20 kpc to the north of 
another elliptical. 

An intriguing case is represented by J1109+3831 whose 
hosting galaxy is a spiral that is located at a projected 
distance of about 20 kpc from an elliptical. A possible 
identification error between optical and radio images has 
been excluded by the analysis of the optical spectrum of 
the companion, which lacks the typical lines displayed by 
active galaxies. 

Among the HFPs identified with galaxies, 3 objects 
(J1058+3353, J1530+2705, and J1602+2646) have an 
optical spectrum in the SDSS DR7. For J1602+2646, 
the optical spectrum seems to be well fitted by a QSO 
template, since a large fraction of the light comes from the 
nuclear region. However, both Fig. [1] and the analysis of the 
diagnostics 0[n]/H^ - OpIIJ/H^ clearly indicate that this 
source is hosted in a galaxy. 

The emission lines detected in these objects (Table [2]) 
are those typical for radio galaxies, showing [O H] A3727, 
[OHI] AA4959, 5007, H^/NH, Yip lines, and [OI] A6300 
in J1058+3353 and in J1602+2646. The [On]/[OIII] line 
ratios indicate low ionization ^ as also found in the optical 
spectra of faint GPS sources (jSnellen et al.lll999l '). 
In the spectrum of each galaxy, the absorption lines associ- 
ated with stellar populations (4300-A G band, 5175-AMg, 
and 5900-ANa D) are clearly visible (Fig.[2l), in particular in 
J1530-I-2705, where the relatively small redshift allows the 



detection of the Call lines. The spectra of the companion 
galaxies (not shown here) have absorption features similar 
to those found in the HFPs, but without prominent [OH], 
[OHI], and H^ emission lines, as expected in non-active 
objects. 

To check whether the galaxies of t he faint HFP sampl e 
follow the Hubble relation found by [Snellen et all (| 199611 ■ 
we added the 12 galaxies with spectroscopic and photo- 
metric redsh i ft in the GPS R-band Hubble diagram of 
'Snel len et al] (|2002l ). The Hubble diagram, where the HFP 
galaxies (squa res) have bee n adde d to the GPS galaxies 
(crosses) from ISnellen et all (|2002l '). is presented in Fig. O 
and it indicates that HFP galaxies have a tight distribution 
in the apparent magnitude-redshift relati on, as also found 
in GPS galaxies (|Snellen et al.lll996l . l2002l ). 



4 RESULTS 

4.1 Spectral properties 

Simultaneous multi-frequency observations carried out at 
different epochs are necessary to monitor the spectral be- 
haviour and variability of high-frequency peaking radio 
sources. Variations in the spectral properties, like peak fre- 
quency, spectral shape and flux densit y, are strong indica- 
tors of the true nature of the source (|Orienti et al.l l2007l : 
iTingav fc de Kool 1120031 '). In young radio sources the spec- 
tral properties are not expected to change, while spectral 
variability is a typical characteristic of beamed radio sources. 
To determine a possible variation of the spectral peak, for 
each source we fitted the simultaneous radio spectrum at 
each epoch with a pure analytical function: 

LogS = a + Log!/ • (b + cLogi^) 

where S is the flux density, v the frequency, and a, b 
and c are numeric parameters without any direct physical 
meaning. We prefer to adopt this f unction instead of 
that used by IStanghellini et aTl (|2009l ) because it better 
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Table 2. Spectral lines of the HFP galaxies with available optical spectrum from the SDSS DR7. Column 1: source name (J2000); Col. 
2: line; Col. 3: line frequency in the observer's frame; Col. 4: line flux density in the observer's frame; Col. 5: equivalent width in the rest 
frame. 



Source 


Line 


A 


'~^linG,obs 

10 " erg s ^ cm ^ 


EWrost 

A 


J1058+3353 


fO 11 


7915.6±1.0 


2.1±0.9 


5.4±0.7 




fO III 


4713. 7±5. 3 


6.3±1.0 


10.0±0.7 




\0 111] 


6273. 6±0. 8 


0.9±0.3 


3.2±0.4 




[OIII] 


6334.2±0.3 


9.9±1.0 


10.0±0.4 






6149.1±0.9 


2.0±0.4 


5.8±0.5 




Ha 


8304.2±6.2 


34.3±1.0 


26.0±0.8 


J1530+2705 


[OI] 


6507.5±0.6 


O.SitO.l 


0.6±0.1 




[Oil] 


3848.3±4.1 


18.3±3.3 


5.5±0.3 




[OIII] 


5120.8±5.1 


1.1±0.2 


1.6±0.2 




[OIII] 


5171.2±0.1 


6.9±0.7 


2.2±0.1 




Hq: 


6779.5±5.0 


6.8±0.3 


2.3±0.1 


J1602+2646 


[OI] 


8643.8±2.4 


24.7±5.2 


26.4±2.6 




[Oil] 


5111.8±0.4 


239.7±15.8 


115.6±2.9 




[OIII] 


6800.7±0.8 


4.7±0.7 


13.2±0.9 




[OIII] 


6865.6±0.5 


23.5±2.1 


29.6±1.1 






6669.1±0.3 


19.0±1.6 


21.9±0.8 




Ha 


9004.7±22.5 


519.6±10.5 


164.7±3.3 



represents the data, providing more accurate values for the 
peak parameters. The best fits to the spectra are shown 
in Fig. |4l and the derived peak frequencies at the various 
epochs are reported in Table ID Statistical errors derived 
from the fit are not representative of the real uncertainty 
on the estimate of the peak frequency. For this reason we 
prefer to assume a conservative uncertainty on the peak 
frequency of 10%. The position of the spectral peak is 
well constrained when the peak occurs at a frequency well 
sampled by the observations, becoming less accurate when 
the frequency coverage is not as appropriate. For example, 
in the case a source has a spectrum peaking at the edges 
of the frequency coverage (i.e. L or K/Q bands), the fit 
provides parameters that are less constrained than in the 
case of sources with the spectral peak occurring around 
5-10 GHz, where both the optically thin and thick emission 
are properly sampled. Spectra poorly constrained are also 
those lacking observations at some frequencies, as in the 
case of J1044+4328, J1052+3355, and J1547+3518 where 
the lack of data either at 1.4/1.7 or 15.3 GHz precluded a 
reliable determination of their peak frequencies. 
By comparing the distribution of the peak frequency of all 
the sources at the various epochs, we do not find remarkable 
differences; the changes are usually within the uncertainties. 
The median value of the peak frequency of the whole sample 
at each epoch has not changed significantly: Up — 5.8 GHz 
in the 1998-1999 epoch (namely the observations of the 
HFP candida tes that have been used for selecting the faint 
HFP sample: [Stanghellini et"ai]|2009l '): i/p = 6.2 GHz in the 
2003-2004 epoch; and fp = 5.4 GHz in the 2006-2007 epoch. 
Among the sources studied here, a few cases show signif- 
icant variation: in 7 sources (J1053-I-4610, J1058+3353, 
J1203-h4803, J1251+4317, J1258-f2820, J1300-K4352, and 
J1616-I-4632) the peak measured in the most recent epoch 
has substantially shifted towards either higher or lower 



frequencies. This result indicates that these sources are part 
of the blazar population, since in young radio sources the 
shift of the peak towards lower frequencies is not expected 
to be so remarkable at least on this short time scale (see 
Section 5). 

We computed the spectral index both below (Qb) and above 
(oa) the peak frequency, by fitting a straight line in the 
optically-thick and -thin part of the spectrum, respectively, 
following the appr oach by iTorniainen et all 12005) and 
lOrienti et al.l (|2007l ). We considered "fiat" those sources 
with both Qb > —0.5 and Qa < 0.5. In a few sources, 
depending on the peak frequency, we could fit either ah or 
Oa only, in order to avoid the flattening near the peak. We 
find that 12 sources (2 galaxies, 9 quasars, and 1 empty 
field) show a flat spectrum during at least one of the 
observing epochs, implying that they are part of the blazar 
population. The fitted optically thick and thin spectral 
indices are reported in Table |3l 

A case that is worth discussion in detail is the quasar 
J1008+2533. The radio spectrum shown by this source 
during two of the observing runs presented here turned 
out to be a composition of two different spectra: convex 
at frequencies below 8.4 GHz, and inverted at higher 
frequencies (see Fig. [4]). This shape is similar to that 
shown by the bright HFP J0927-I-3902. In J0927-I-3902 
the two-compone nt spectrum is e xplained by its core-jet 
structure (see e.g. lOrienti et al1l2QQ6a i: at frequencies below 
1 GHz the spectrum is dominated by the emission from 
the jet, while at higher frequencies the contribution from 
the self-absorbed core becomes more important, becoming 
the dominant emission above ~ 10 GHz. Such a scenario 
is well supported by pc-scale m orphological informa tion 
by muhi- frequency VLBI data (|Alberdi et al.l bOOd ). A 
similar explanation may apply to the case of the faint HFP 
J 1008-1-2533. Another possibility is that in our new epochs 
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Table 3. Multi-frequency VLA flux density of 57 sources from the faint HFP sample. Column 1: source name (J2000); Col. 2: optical 
identification from the SDSS DR,7: G=galaxy, Q=quasar; EF=empty field; Col. 3: R magnitude; Col. 4: redshift. An "f" indicates a 
photometric redshift from the SDSS DR7; Col. 5: the observing code from Tablc[T] Cols. 6 - 14: flux density at 1.4, 1.7, 4.5, 5.0, 8.1, 8.4, 
15.3, 22.2, and 43.2 GHz respectively; Cols. 15 and 16: the spectral index computed below and above the spectral peak, respectively. 
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Table 3. Continued 



Source ID mag z code 51.4 Si. 7 S4.5 S5.0 Ss.i Ssa S15.3 522.2 ^43. 2 cth 

R mjy mjy mjy mjy mjy mjy mjy mjy mjy 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) (13) (14) (15) (16) 



J1240+2323 


G 


21.0 


0.38f 


a 


27 


28 


52 


53 


59 


60 


54 


49 




-0.6 


0.2 










c 


24 




58 


60 


65 


65 




56 




-0.6 




J1240+2425 


Q 


16.9 


0.831 


a 


78 




60 


58 


45 


44 


32 


27 






0.4 










c 


57 




58 


56 


45 


44 




17 






0.5 


J1241+3844 


Q 


21.5 




c 


20 


21 


23 


25 


20 


20 




13 




-0.2 


0.4 


J1251+4317 
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18.7 


1.453 


a 






54 


57 


80 


82 


112 


116 




-0.4 












c 


25 




72 


80 


138 


142 




105 




-1.0 


0.3 


J1258+2820 


Q 


19.2 




a 


25 


32 


51 


52 


51 


50 


43 


38 




-0.6 


0.2 










c 


23 


34 


45 


45 


54 


54 




54 




-0.3 




J1300+4352 


Q 


19.6 




c 


140 


155 


116 


115 


97 


95 




81 






0.2 


J1309+4047 


Q 


18.9 


2.910 


a 


37 




139 


131 


118 


115 


76 


40 




-1.1 


0.8 










c 


34 


51 


132 


133 


113 


110 




33 




-1.1 


0.9 


J1319+4951 


Q 


19.1 




a 


25 




49 


49 


39 


39 


25 


21 




-0.5 


0.5 










c 


22 


30 


50 


49 


43 


43 




20 




-0.7 


0.6 


J1321+4406 


Q 


21.2 




a 


62 




73 


74 


71 


70 


58 


50 




-0.1 


0.3 










c 


65 


73 


78 


77 


75 


74 




38 




-0.1 


0.4 


J1322+3912 


Q 


17.5 


2.985 


a 


117 




227 


223 


181 


176 


118 


86 




-0.6 


0.6 










c 


116 


135 


200 


196 


150 


146 




56 




-0.5 


0.8 


J1330+5202 


G 


20.7 


0.58f 


c 


91 


103 


176 


180 


180 


177 




156 




-0.4 


0.1 


J1336+4735 


Q 


19.7 




a 






61 


61 


52 


51 


26 


22 






0.7 










d 


26 




67 


67 


57 


56 




30 






0.5 


J1352+3603 


G 


18.0 


0.27f 


a 


65 


70 


102 


103 


95 


93 


65 


37 




-0.4 


0.7 










d 


59 




97 


97 


84 


82 




34 




-0.4 


0.7 


J1420+2704 


Q 


20.3 




a 


14 




55 


57 


55 


53 


34 


25 




-1.1 


0.6 










d 


14 




62 


63 


58 


56 




26 




-1.2 


0.6 


J1436+4820 


EF 






a 


20 




72 


72 


61 


58 


32 


13 




-1.0 


1.1 


J1459+3337 


Q 


16.6 


0.645 


a 


22 


30 


195 


221 


403 


415 


470 


435 




-1.3 




J1530+2705 


G 


14.1 


0.033 


a 


13 


14 


43 


46 


45 


43 


17 


13 




-1.0 


0.8 










d 


11 


12 


48 


49 


42 


41 




23 




-1.1 


0.5 


J1547+3518 


Q 


21.2 




a 






53 


57 


67 


67 


74 


72 




-0.3 












d 


12 




49 


51 


57 


58 




74 




-0.6 




J1602+2647 


G 


18.4 


0.372 


a 


39 




131 


148 


238 


244 


259 


224 




-0.8 




J1613+4223 


Q 


19.8 




b 


42 




206 


201 


119 


110 


39 


14 




-1.4 


1.7 










d 


40 




201 


194 


114 


105 




13 




-1.4 


1.7 


J1616+4632 


Q 


19.3 


0.950 


d 


88 




141 


144 


148 


147 




129 




-0.3 


0.1 


J1617+3801 


Q 


19.0 


1.607 


a 


23 




71 


75 


99 


100 


95 


67 




-0.8 


0.4 


J1624+2748 


EF 






a 


20 




109 


118 


175 


177 


198 


173 




-1.0 




J1651+3417 


EF 






b 






46 


48 


58 


58 


49 


22 


16 




0.9 


J1702+2643 


Q 


17.2 




b 


35 


39 


41 


42 


49 


50 


61 


65 


73 


-0.2 




J1719+4804 


Q 


15.3 


1.084 


a 


75 




126 


135 


158 


157 


115 


85 




-0.4 


0.6 










b 


66 




122 


127 


150 


148 


113 


73 


30 


-0.5 


0.7 










d 


60 




112 


114 


100 


97 




43 




-0.5 


0.6 



a flare from a self-absorbed knot in the jet occurred, causing 
an increase of the flux density at high frequencies, as also 
found in blazar objects. The lack of pc-scale morphology 
information does not allow us to unambiguously determine 
the origin of this complex spectrum. 



4.2 Variability index 

Following the ap proach 



i^ollowmg tne ap proacn by iTinti et al.1 lj2005f ) and 
lOrienti et aP (|2007h we investigate the presence of flux den- 
sity variability, by computing the variability index V: 



V = 



_j_ \ ^ {Si — Sif 



m ^—1 

2 = 1 



(1) 



where Si is the flux density at the j-th frequency measured 
at one epoch, Si is the mean value of the Hux density 
computed averaging the flux density at the i-th frequency 
measured at all the available epochs, Oi is the rms on 
5"; — Si, and m is the number of sampled frequencies. We 
prefer to compute the variability index for each new epoch 
(see Table [4| instead of considering all the epochs together 
in or der to bette r dete ct the presence of a possible burst. 
As in iTinti et al.l (|2005l ). we consider variable those sources 
with a variability index V > %. 



8 M. Orienti et al. 



V 



J0736+4744 Q 




V (GHz) 

J0821+3107 Q 



1000.0 0.1 



V (GHz) 



V 


J0955+3335 Q 








&k \\ 
'I \^ 






V (GHz) 



J1020+2910 EF 




H 



J0754+3033 Q 




H 



J0804+5431 G 



H 



JOS 19+3823 Q 




V (GHz) 



V (GHz) 



V (GHz) 



100 1000,0 



H 


J0905+3742 EF 




H 


J0943+5113 G 




H 


J095 1+3451 G 






100 




'A 






/ X 
/ \ 
















* i 






10 





























V (GHz) 



V (GHz) 



V (GHz) 



H 


J 1002+5701 EF 




V 


J1004+4328 EF 




1 I 


;mJy) 








jit ^'A 












10 







J1008+2533 Q 




V (GHz) 

J 1020+4320 Q 




V (GHz) 

J1025+2541 G 




H 



V (GHz) 

J 1035+4230 Q 




V (GHz) 



V (GHz) 



V (GHz) 



V (GHz) 



Figure 4. Radio spectra of the 57 candidate HFPs from the "faint" HFP sample observed with the VLA during the observ ing runs 
presented in this paper. Asterisks and a dash-dotted hne refer to the first epoch observations (1998-2000. IStangheUini et al.|[200gt) : crosses 
and a dashed line refer to epoch a (2003); diamonds and a solid line refer to epoch b (2004); squares and a dotted line refer to epochs 
c,d (2006-2007). 
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Figure 4. Continued. 
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Figure 4. Continued. 
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Figure 4. Continued. 



From the comparison of the 
properties and variability we find: 



multi-epoch spectral 



- 24 objects maintain the convex spectrum without showing 
significant flux density variability {V < 9). They are 
labelled "H" in Column 9 of Table H 

- 21 sources preserve the convex spectrum at the various 
epochs, although with some amount of fiux density variabil- 
ity (1/ > 9). They are labelled "V" in column 9 of Table H 

- 12 sources with a co nvex spectrum in the first epoch 
jStanghelUni et al.l [2OO9I ') show a flat spectrum at least in 
one of the subsequent epochs. They are marked with "F" 
in Column 9 of Table g] 



4.3 Radio properties versus optical identification 

The optical identification of the radi o sources in sam- 
ples of CSS and GPS objects (e.g [Fanti et all Il990l: 
Stanghellini et al.l 1 19981 : ISnellen et all 1 19981 : iFanti et all 
2001f ) showed that in "faint" samples there is a higher 



fraction of objects identified with galaxies with respect 
to "bright" samples, and the fraction of galaxies seems 
to be anti-correlated with the peak frequency, i.e. in 
CSS samples a higher percentage of radio sources are 
hosted in galaxies than in GPS and HFP samples. This 
is easily seen in Fig. [SI where we have plotted the peak 
frequency versus the peak flux density for the source s 
from the bright HFP s ample JPallacasa et al.l l200Cll '). 
the faint HF P sample ll Stanghellini et al.l 120091 ') . the 
GPS sa mple (IStangheUini^ all 1 19981 ). "the faint GPS 
sample (| Snellen et al. 1998h : tlie half-Jansky GPS sam- 
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Figure 2. Optical spectra from the SDSS DR7 of tlie galaxies 
J1058+3353 (top), J1530+2705 {center), J1602+2642 (bottom). 



Figure 3. The Hubble diagram of HFP galaxies (squares) and 
GPS galaxies (crosses, Snellen et al. 2002). The dotted line 
represents th e Hub ble relation for GPS galaxies as found by 
[Snellen et al] jiggd) . 



pie (ISnellen etj] |2002| ). and the B3-VLA CSS sample 
(Fanti et al. 200 ll). In the latter sample only objects with 
spectra peaking above 100 MHz have been considered, since 
for sources with a peak at lower frequency, the turnover 
frequency could not be reliably constrained. In this plot, 
radio sources identified with galaxies (squares) mainly 
occupy both the top left part of the panel, i.e. low peak 
frequency and high peak flux density, and the bottom right 
panel, i.e. high peak frequency (usually below 10 GHz) 
but low peak flux density. On the other hand, quasars 
(circles) are mostly found in the left part of the plot, with 
peak frequencies above a few GHz, and peak flux densities 
that span almost three orders of magnitude. Radio sources 
without an optical identiflcation (triangles) are found in 
the same regions as galaxies, suggesting that the majority 
of these objects share the same properties of the galaxies, 
but they are fainter probably because they are at higher 
redshift. It is worth noting that the criteria that have been 
used to date to select radio source catalogues cannot pick 
up objects with both low peak frequency and low peak flux 
density (bottom left panel), indicating that we miss the 
population of faint CSS objects. 

Statistical studies on the correlation between the radio 
characteristics and the optical counterpart indicate that 
those objects hosted in a galaxy have the typical properties 
of young radio sources (i.e. symmetric radio structure 
and no spectral variability), whereas those identifled with 
quasars are more similar to flat-spectrum radio objects. 
This suggests that there is a dichotomy between the optical 
identiflcation and the radio properties: quasars are more 
likely part of the flat-spectrum blazar population, while 
galaxies are likely associated with genuinely young radio 
sources. In the case of the sources from the faint HFP 
sample studied in this paper, the comparison between their 
optical identiflcation and the radio variability showed that 
the majority of quasars have a variable spectrum, while 
a smaller fraction maintain the convex spectrum without 
variability (33% H, 39% V, and 28% F). In the case of 
galaxies, we found that the fraction of variable objects is 
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still larger than those without spectral variability (34% 
H, 50% V, and 16% F). For the sources still lacking an 
optical identification we found that the majority do not 
show significant spectral variability (75% H, 17% V, 8% F). 
Among the 12 sources with a fiat radio spectrum, we 
found that two objects (J1058+3353 and J1330+5202) are 
identified with a galaxy, and one (J1020+2910) lacks an 
optical identification. When we compare the variability 
properties of sources with different optical identification by 
means of the Student's t-statistic we find that there is a 
significant difference (>90%) between galaxies and quasars, 
indicating that the majority of galaxies and quasars are 
associated with different radio source populations. 



5 DISCUSSION 

The anti-correlation found betwe en the projected lin- 
ear size and the p e ak frequency l O'Dea fc Baum I Il997l : 
iBicknell et"all 1 19971 : ISnellen et all \200&i implies that the 
sources with the spectral peak occurring above a few 
GHz should represent the population of the smallest 
radio sources whose radio emission has recently turned 
on. Samples of high-frequency peaking objects have been 
selected by choosin g sources with an inverted radio spec- 
trum up to 5 GHz (|Dallacasa et aLlliood : iTorniainen et all 
I2OO5I : IStanghellini et al.l l2009l ). i.e. the highest observing 
frequency where a large area survey is presently available. 
However, due to the selection criteria these samples have 
been found to comprise both young radio sources and 
flaring flat-spectrum objects selected during particular 
phases of their spectral variability, for example when their 
radio emission is dominated by a knot in the jet. The 
study of flux density and spectral variability based on 
repeated simultaneous multi-frequency observations has 
proved to be an ideal tool in discriminating the different 
nature of the sources. It was found that in samples of bright 
objects, where there is a high incidence of sources optically 
identifled with quasars, fiat-spectru m blazar objects repre- 
sent the dominant population (e.g. Torniainen et al.1 l2007l : 
lOrienti et"al]l2007l : Ijauncev et al.ll2003l l. A higher incidence 
of genuinely young radio sources is expected in samples of 
faint objects where the majority of radio sources should be 
hosted in galaxies and boosting effects are supposed to play 
a minor role. 

The optical identification of the sources studied in this 
paper by means of the SDSS DR7 indicates that 21% of 
objects are hosted in galaxies, i.e. similar to the fraction of 
galaxies in the bright HFP sample. However, in the faint 
sample, another 21% of objects lack optical identification 
and thus a reliable comparison between the two samples 
cannot be done. The analysis of the optical images of 
the galaxies hosting HFP pointed out the presence of 
companions around 6 HFP candidates, indicating that 
young radio galaxies, like powerful exten ded radio sources, 
are in groups, as previously suggested bv lStanghellini et al.l 
(| 19931 ). indicating a continuity between compact young 
objects and the population of classical radio galaxies 
(jO'Deaet al.lll996l ). 

Although in 4 galaxies the presence of companions is 
suggested by photometric information only, in J1530-f 2705 



and J 1602-1-2646 the association is made by spectroscopic 
redshift. The companion galaxies are located within a pro- 
jected distance of about 150 - 200 kpc from the target which 
usually is hosted in the brightest elliptical at the group 
centre. A peculiar case is represented by J1109-I-3831, whose 
parent galaxy seems to be a spiral that is interacting with 
an elliptical. Young radio sources are normally associated 
with ellipticals. The case represented by J1109-I-3831 may 
be explained by the possible interaction between the hosting 
spiral and the companion that may have triggered the radio 
emission. The small redshift of J 1530-1-2705 enabled us to 
identify the morphology of its brightest companions that 
turned out to be barred spirals. This group resemble s that 
of the bright HFP J0655-f4100 (|Orienti et al.ll2006bl ). and 
in both cases the HFP is hosted by the central elliptical 
galaxy at the group centre. The presence of companion 
galaxies in the environment of galaxies hosting young 
radio sources suggests that the onset of the radio emission 
may be triggered by merger or interaction events that 
occurred not long ago. This scenario is supported by 
the proximity of the companions in J0804+5431 and in 
J1109-I-3831, although observations to establish a physi- 
cal interaction are needed to unambiguously verify this idea. 

From the analysis of the multi-epoch radio spectra of 
the sources in the faint HFP sample, we find a high frac- 
tion of objects displaying some level of variability. This 
result does not imply that all these sources are part of 
the blazar population. In fact, changes in the radio spec- 
trum may be a direct consequenc e of the source expan- 
sion (e.g. iTingav fc de Kool I I2OO3I ). In newly born radio 
sources, the evolution timescales can be of the order of a 
few tens of years. Changes in the radio spectrum of such 
young objects can be appreciable after the short time (5-8 
years) elapsed between the first and last observing run. A 
clear example is represented by the faint HFP J1459-I-3337 
(|Orienti fc Dallacasa IbOOSd ). This HFP showed a steadily 
increasing flux density at 1.4 and 5 GHz, in the optically 
thick regime, and its spectral peak s hifted f rom 30 GHz 
down to 12 GHz in about 10 years (jEdge et al.lil996 ). This 
behaviour is consistent with the flux density and spectral 
evolution of a young object, with an age of about 50 years, 
undergoing adiabatic expansion. 

In the presence of adiabatic expansion of a homogeneous 
synchrotron source, the radio spectrum undergoes a shift 
towards low frequencies. In the optically-thick regime this 
means that at a given frequency the flux density S increases 
with time (see Pacholczyk 1970, Orienti et al. 2007, and Ori- 
enti et al. 2008b for a detail analysis of the radio spectrum 
evolution) : 



5*1 — So 



to + At 
to 



(2) 



where So and 5*1 are the flux densities at the time to and 
to + At, respectively. On the other hand, the spectral peak 
i/p moves to lower frequencies: 



to 



to+At 



(3) 



where Upfi and Vp^i are the peak frequency at the time to 
and to + At respectively. 
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Table 4. Peak frequency and flux density variability. Column 1; source name (J2000); Col. 2: peaJc frequency in GHz of the first epoch 
(1998-2000, Stanghellini et al. 2009); Cols 3, 4, and 5: peak frequency in GHz of epoch a (2003), b (2004), and c, d (2006-2007); Cols. 6, 
7, and 8: variability index V computed for the epoch a, b, c,d respectively; Col. 9: the classification of the source spectrum (V=variable, 
H=genuine HFP, F=fiat; see Section 3). An asterisk indicates that the peak frequency is not reliable due to poor frequency coverage. 
for the source J1008+2533 we report the peaJs of the lowest part of the spectrum. 
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Table 4. Continued. 



Source i^opi i'cp2 I'opS i'cp4 Va Vfc Vc,d Var. 
(1) (2) (3) (4) (5) (6) (7) (8) (9) 
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In the optically-thin regime, the flux density at a given fre- 
quency decreases with time: 



5*1 = 5*0 



(4) 



where S is the spectral index of the electron energy dis- 
tribution N{E) (X that originates the radio emission. 
It is clear from this relationship that substantial variation 
in the optically-thick regime can be revealed in case that 
At is a non-negligible fraction of the total source age (to)- 
In our case, At ~ 8 years and it can produce detectable 
flux density variation for sources with to ~ 100 years. 
Among the sources with some changes in the spectral 
properties compatible with such a scenario, other 7 sources 
(J0754-h3033, J0955-H3335, J1004-H4328, J1025-H2541, 
J1052-F3355, J1322-F3912, and J1547-H3518) in addition to 
J1459-I-3337 show a spectral and flux density variability that 
may be explained in terms of source expansion. Although 
in J1459-I-3337 there are several indicators s upporting this 
interpretation (lOrienti &: Dallacasa I l2008d ). in the case 
of the other sources mentioned here this assumption is 
based on the flux density and peak variation on a small 
time range only. Additional observations spanning a longer 
time interval, together with information on the pc-scale 
morphology are necessary in order to reliably constrain the 
source nature. 
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Figure 5. The peak frequency versus the peak flux density for 
the radio sources from the B3- VLA CSS sample jFanti et al.l 
l200j) . the faint GPS sample (Snellen ct al.' 19981), the half- 
Jansky GPS samp le (Snellen et al. 2002), the GPS s ample from 
IStanghellini et al.l (|l998), the bright (Dallacas a et ai]|2000! ) and 
the faint HFP samples ( Stan glicllin i e t al . 200!^ ). Circles, squares 
and triangles refer to quasars, galaxies and empty field, respec- 
tively. 



6 CONCLUSIONS 

We presented simultaneous multi-frequency VLA obser- 
vations of 57 sources from the faint HFP sample, carried 
out at various epochs. From the comparison of the spectral 
properties we found that 24 objects (4 galaxies, 11 quasars, 
and 9 empty fields) preserve their convex spectrum without 
showing any evidence of flux density variability. Of the 
remaining sources, 12 objects (2 galaxies, 9 quasars, and 1 
empty field), selected on the basis of the i r conv ex spectrum 
in the first epoch by IStanghelhni et all jlooi), turned out 
to show a flat spectrum in one of the subsequent observing 
epochs. The remaining 21 sources (6 galaxies, 13 quasars, 
and 2 empty fields) possess high levels of variability, 
although still displaying a convex spectrum. However, 
among these variable sources we found that in 8 objects the 
changes in their spectra are consistent with what expected 



if the source is undergoing adiabatically expansion. This 
implies that out of the 57 sources studied in this paper, 32 
objects (56%) can still be considered young radio source 
candidates. The remaining 25 sources (44%) are part of 
the flat-spectrum blazar population, indicating that also 
in samples of faint radio sources, where boosted effects are 
thought to play a minor role, a large fraction of sources are 
represented by flaring objects. 

The analysis of the optical images of the HFPs hosted by 
galaxies pointed out the presence of companion galaxies 
in the target environment, supporting the idea that young 
radio sources reside in groups. The parent galaxy is usually 
the brightest elliptical at the group centre with the ex- 
ception of two sources. In J0804-I-5431 the galaxy hosting 
the HFP is at the periphery of the group, and it seems 
interacting with a close elliptical. A surprising result is 
represented by the HFP J1109+3831 that is hosted in a 
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spiral that seems to be interacting with a close elliptical. 
The fact that young radio sources reside in groups support 
the idea that the interactions occurring between the galaxies 
are at the origin of the radio emission. 
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